Dori-Bachash et al. Genome Biology 2012, 13:R1 14 
http://genomebiology.com/201 2/1 3/1 2/R1 1 4 



^ Genome Biology 



RESEARCH 



Open Access 



Widespread promoter-mediated coordination of 
transcription and mRNA degradation 



Mally Dori-Bachash 1+ , Ophir 5halem 1+ , Yair S Manor 1 , Yitzhak Pilpel 1 * and Itay Tirosh 



1,2* 



Abstract 

Background: Previous work showed that mRNA degradation is coordinated with transcription in yeast, and in 
several genes the control of mRNA degradation was linked to promoter elements through two different 
mechanisms. Here we show at the genomic scale that the coordination of transcription and mRNA degradation is 
promoter-dependent in yeast and is also observed in humans. 

Results: We first demonstrate that swapping upstream c/s-regulatory sequences between two yeast species affects 
both transcription and mRNA degradation and suggest that while some c/s-regulatory elements control either 
transcription or degradation, multiple other elements enhance both processes. Second, we show that adjacent 
yeast genes that share a promoter (through divergent orientation) have increased similarity in their patterns of 
mRNA degradation, providing independent evidence for the promoter-mediated coupling of transcription to 
mRNA degradation. Finally, analysis of the differences in mRNA degradation rates between mammalian cell types 
or mammalian species suggests a similar coordination between transcription and mRNA degradation in humans. 

Conclusions: Our results extend previous studies and suggest a pervasive promoter-mediated coordination 
between transcription and mRNA degradation in yeast. The diverse genes and regulatory elements associated with 
this coordination suggest that it is generated by a global mechanism of gene regulation and modulated by gene- 
specific mechanisms. The observation of a similar coupling in mammals raises the possibility that coupling of 
transcription and mRNA degradation may reflect an evolutionarily conserved phenomenon in gene regulation. 



Background 

Gene expression patterns are regulated by a multitude of 
processes, and while these processes have traditionally 
been studied in isolation, numerous instances of cross-talk 
between them have been revealed in recent years [1-5]. 
These cross-talk events influence the efficiency and 
dynamics of gene regulation and reflect the high inter- 
connectivity between regulators of distinct processes. 
While initial studies focused on the cross-talk between 
successive processes that occur in the same cellular com- 
partment (for example, co-transcriptional mRNA proces- 
sing), recent evidence points to coupling between 
transcriptional events in the nucleus and mRNA regula- 
tion in the cytoplasm. For example, Rpb4/7 binds to 
mRNAs co-transcriptionally and shuttles with the mRNA 
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to the cytoplasm where it influences mRNA degradation 
and translation [6-8]. Large-scale analyses in yeast have 
demonstrated a global coordination between transcription 
and mRNA degradation in response to environmental 
[9,10] and genetic [11] perturbations and during evolution 
[12], suggesting that the mechanistic coupling through 
Rpb4/7 and other components has a global impact on 
gene regulation. 

Coupling of transcriptional induction with enhanced 
mRNA degradation generates a peaked response that is an 
important feature of dynamic biological systems that must 
respond quickly but transiently as they adapt to environ- 
mental changes [10,13,14]. However, we recently demon- 
strated, through comparison of yeast species grown in rich 
media [12], that this coupling is also evident in the steady- 
state regulation of gene expression. Thus, it is possible 
that the transcription-degradation coupling has evolved in 
order to facilitate peaked responses but subsequently 
became a fundamental aspect of regulation. The presence 
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of this global coupling in other species (besides yeast) 
remains unclear. 

Evolutionary changes in transcription and mRNA 
degradation were coupled both through trans and 
through cis mutations, as determined from analysis of 
interspecific hybrids [12]. Coupling through trans muta- 
tions presumably involves changes in the activity of pro- 
tein complexes that regulate both transcription and 
degradation, such as Rpb4/7 and Ccr4-Not [15]. Cou- 
pling by ds-mutations was observed at more than a hun- 
dred genes, but the mechanism of ds-coupling mutations 
is unknown. We speculated that promoter cis mutations 
could affect the recruitment of transcription complexes 
that not only regulate transcription but also prime the 
newly transcribed mRNA for future regulation in the 
cytoplasm [16]. Interestingly, two recent studies demon- 
strated promoter-mediated effects on mRNA degradation 
through binding of Rapl or through cell-cycle regulators 
[17,18]. Yet these specific mechanisms are unlikely to 
account for our previous observations [12] since the cis- 
coupled genes we discovered were underrepresented with 
Rapl targets and with cell-cycle periodic genes (not 
shown). This suggests the existence of additional 
unknown coupling mechanism(s) through c«-regulatory 
sequences, which may involve promoter, UTR or open 
reading frame elements. 

Here, we validate our previous prediction of cw-depen- 
dent coupling by swapping upstream czs-regulatory 
sequences (promoters and 5' UTRs) between Saccharo- 
myces cerevisiae and Saccharomyces paradoxus and mea- 
suring the effect of swapping on mRNA levels and mRNA 
degradation. We further demonstrate the role of promoters 
in controlling mRNA degradation by showing that yeast 
genes that share a promoter tend to have similar mRNA 
degradation profiles. Based on these results we propose 
that a general promoter-mediated mechanism couples 
transcription to mRNA degradation of a large fraction of 
yeast genes and that this mechanism may be modulated by 
various cis and trans elements. Finally, we present evidence 
for a similar coupling in mammalians through analysis of 
differences in mRNA degradation between different cell 
types or species. 

Results 

Swapping orthologous sequences reproduces 
interspecies expression divergence and 
transcription-degradation coupling 

In previous work, we used interspecific hybrids to predict 
cw-dependent changes in gene expression for hundreds 
of genes, and a ris-dependent coupling of transcription 
and mRNA degradation. To validate these predictions 
and to examine whether the relevant czs-acting mutations 
may reside at the promoter or 5' UTRs, we generated 
S. cerevisiae strains in which upstream czs-regulatory 



sequences (200 to 700 bp of the sequence immediately 
upstream to the start codon of the gene of interest) 
have been replaced by the orthologous sequence from 
S. paradoxus, which are approximately 80% identical 
(Figure la; Table SI in Additional file 1). For most genes 
examined (18 of 34), swapping affected mRNA levels and 
reproduced the previously predicted interspecies cis- 
dependent differences (Figure lb). 

To test whether the differences in mRNA levels reflect 
differences in transcription, mRNA degradation or both, 
we focused on nine of the genes examined above that 
were previously identified as having a czs-dependent cou- 
pling of transcription and mRNA degradation. We 
arrested transcription (with 1,10-phenanthroline) of 
strains with the wild type (WT) and with the swapped 
(Swapped) regulatory sequences of these genes, and mea- 
sured the mRNA log 2 -ratio (Swapped/WT) at four time 
points (Figure lc). The log 2 -ratio at the zero time point 
(before transcriptional arrest) reflects the difference in 
steady-state mRNA levels, while the change in that ratio 
along the time course reflects differences in mRNA 
degradation rates: if the swapped sequences increase 
mRNA degradation rates, then the mRNA log 2 -ratio 
(Swapped/WT) should decrease linearly across the time 
course (marked in red at Figure lc), while if mRNA 
degradation decreases by the swapping then the mRNA 
log 2 -ratio should increase along the time course (marked 
in green). Strikingly, swapping affected both the steady- 
state mRNA levels and the rates of mRNA degradation 
for most genes examined (7 of 9), and in all of these 
cases the change in mRNA levels was opposite to that 
expected from the change in mRNA degradation: 
increased degradation was associated with increased 
(rather than decreased) steady state mRNA levels, and 
vice versa. This indicates that changes in mRNA degrada- 
tion were accompanied by changes in transcription that 
had an opposite effect on mRNA levels, and that the 
changes in transcription were larger and thus determined 
the direction of change in mRNA levels. Four of these 
genes were chosen for validation (see Materials and 
methods), and their results were reproduced in additional 
experiments with eight time points (Figure Id). 

As a control, we examined the effect of the selection 
marker that is inserted upstream of the swapped c/s-regu- 
latory elements (Figure la). We found that inserting 
the markers at the same location, but without swapping 
the as-regulatory sequences, had moderate effects on the 
expression levels of the regulated genes that did not 
account for most changes observed in the swapped strains 
(Figure Sla in Additional file 1). Assuming that the 
inserted marker and the swapped cis-regulatory region 
exert independent effects on the regulated genes, we sub- 
tracted the estimated marker effects from the observed 
effects of 17 swap strains described above. This analysis 
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Figure 1 Swapping of upstream regulatory sequences affects transcription and mRNA degradation, (a) Scheme describing the generation of 
5. cerevisiae Swap strains by replacing the native S. cerevisiae promoter and 5' UTR with the orthologous sequences from S. paradoxus. Note that for 
each strain we modulate the regulatory region of only one gene while all other loci are identical to the wild-type (WT) S. cerevisiae (b) Log 2 -ratio of 
mRNA levels in the swapped versus WT strains are shown for 34 genes with a predicted c/s-dependent interspecies expression difference. Colors 
represent the concordance between the swapping effect and the inter-species expression difference. Expression changes of less than 40% were 
considered not significant, as the interspecies expression differences were approximately two-fold, (c) Log 2 -ratio of mRNA levels in the Swapped 
versus WT strains after transcriptional arrest. Colors indicate whether mRNA degradation rates are increased (red), decreased (green) or not significantly 
affected (black, P > 0.05). The average values from three technical replicates are shown along with regression lines. Error bars show standard deviation 
among three technical replicates, while both values are shown for experiments with two biological replicates (AAT2, MKC7, OXP1 and YML108W). 
(d) Same as (c) for a validation experiment of four selected genes with eight time points. 
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slightly increased the proportion of genes in which the 
swapped region reproduced the predicted as-dependent 
interspecies differences (Figure Sla in Additional file 1). 
Surprisingly, for the two genes that were most sensitive to 
the selection marker (Pex32 and Gal83) the swapped 
strains also had an effect on mRNA degradation, suggest- 
ing that these effects could be due to the selection marker 
instead of the swapped regulatory region. Indeed, analysis 
of control strains (with the marker but without swapping 
the regulatory sequences) indicated that, for some genes, 
the selection marker affected both transcription and 
mRNA degradation (Figure Sib in Additional file 1). Thus, 
although the mechanism by which insertion of the selec- 
tion marker affected the expression of these genes is 
unknown, this mechanism also coupled the transcription 
and mRNA degradation of some genes. A possible expla- 
nation is that these genes, which were chosen based on 
previous evidence of transcription-degradation coupling, 
are regulated by a mechanism that links their transcrip- 
tional activity to mRNA degradation and therefore that 
various means of modulating their transcriptional activity 
could induce an effect on mRNA degradation. 

Detailed analysis reveals coupled and uncoupled effects 
on transcription and mRNA degradation 

Two genes were selected for a more detailed analysis (see 
Materials and methods). For each of the two genes, we gen- 
erated a series of sequential S. cerevisiae strains where 
different fractions of the respective upstream regulatory 
sequence have been swapped by the orthologous region 
from S. paradoxus (Figure 2a). Each pair of 'adjacent' strains 
therefore differs by the species-of-origin (S. cerevisiae versus 
S. paradoxus) of a relatively small sequence element, and 
comparison of these strains uncovers the effect of mutations 
in this sequence only. 

We first analyzed six swap strains of MRU and found 
that they are divided into two clusters that differ from the 
control strain (Figure 2b,c). The first cluster had a large 
decrease in mRNA levels (>2-fold) and a more subtle 
decrease in mRNA degradation. This cluster includes 
three strains with swapping of 107, 168 and 197 bp (out of 
the complete 300 bp swapped segment), indicating that 
the causal sequence is in the 107 bp upstream-most por- 
tion of the swapped segment. Consistent with previous 
studies [19,20], 5' RACE analysis (Figure S2 in Additional 
file 1) mapped the transcription start site of MRU to 
approximately position 250 within the swapped sequence 
(that is, 5' UTR of approximately 50 bp or even less) in the 
WT and Swapped strains, indicating that the causal 
sequence is a promoter element that couples a transcrip- 
tional effect with a (weaker) mRNA degradation effect. 

In contrast, the second cluster showed an additional 
decrease of approximately two-fold in mRNA degradation 
and a more subtle increase in mRNA levels (compared to 



the first cluster). The increase in mRNA levels is expected 
from the decreased rates of mRNA degradation, suggest- 
ing that mRNA degradation is the primary effect and that 
transcription was either unaffected or only mildly affected 
by the additional swapped sequences. This cluster includes 
two strains with swapped segments of 225 and 300 bp, 
indicating that the causal sequences that distinguish this 
cluster from the first cluster are in the 197 to 225 bp 
region. This region is also upstream of the mapped tran- 
scription start site) and therefore reflects a promoter ele- 
ment that controls mRNA degradation. 

Next we analyzed the swap strains of YLR326W (Fig- 
ure 2d,e; Figure S3 in Additional file 1). Swapping of the 
entire 470 bp upstream region increased mRNA levels 
by approximately four-fold and mRNA degradation by 
approximately three-fold. The series of swap strains 
indicated a large effect of an upstream region (positions 
200 to 300) and smaller effects of multiple other regions 
between positions 300 to 471 (Figure 2e). 5' RACE ana- 
lysis (Figure S2 in Additional file 1) mapped the tran- 
scription start site to two positions corresponding 
approximately to positions 340 and 400, indicating that 
the upstream region with the large effect is a promoter 
element. This promoter element appears to increa- 
semRNA degradation by approximately two-fold and 
transcription by approximately four-fold, thereby gener- 
ating a net increase of two-fold in mRNA levels. The 
rest of the swapping effect, which was not accounted by 
this promoter element, was spread among eight distinct 
segments (positions 300 to 471) that collectively account 
for a 2-fold increase of mRNA level and a 1.5-fold 
increase of mRNA degradation. Each of these elements 
had effects that may be too weak to interpret individu- 
ally (up to 30% changes in mRNA levels and degrada- 
tion by each segment) as our measurements could have 
been affected by technical and biological variability. 
Nonetheless, we note that at least one of these segments 
(positions 386 to 398) appears to affect both transcription 
and mRNA degradation (each by approximately 30%) as 
these effect were statistically significant (P < 0.05; see 
Materials and methods) and were consistent with two 
independent experiments in which we swapped the con- 
taining segment 360 to 414 and observed a similar effect 
(Figure S3 in Additional file 1). This segment is found 
between the two distinct transcription start site positions 
of YLR326W and thus could control transcription and 
mRNA degradation either as a promoter or as a 5' UTR 
element. 

Gene-pairs that share a promoter have similar patterns 
of mRNA degradation 

To provide further statistical genome-wide support to the 
emerging notion that promoters might regulate mRNA 
degradation, we turned to examine the prediction that 
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Figure 2 Analysis of smaller swapped segments reveals multiple distinct effects on transcription, mRNA degradation and their 
coupling, (a) Scheme describing the generation of a series of sequential 5. cerevisiae Swap strains by replacing the different fractions of the 
native S. cerevisiae promoter and 5' UTR with the orthologous sequences from 5. paradoxus. (b,d) Log 2 -ratio of mRNA levels in the Swapped 
versus WT strains after transcriptional arrest of MRU (b) and YLR326W (d). Error bars reflect standard error among two or four biological replicates 
(see Materials and methods). (c,e) Differences in mRNA degradation rates estimated using the linear least square fits shown in (b,d). Swap strains 
were classified into clusters by starting from a single cluster for all strains and separating consecutive strains into distinct clusters if they are 
significantly different (P < 0.05) in mRNA levels and/or mRNA degradation; the significant changes are noted by arrows. 
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pairs of genes that are co-regulated by the same promo- 
ter elements will have similar patterns of mRNA degrada- 
tion. Due to the compactness of the yeast genome the 
intergenic regions are often very short (typically lower 
than 1 kb) and hence neighboring gene pairs, with diver- 
gent orientation [21] (Figure 3a), might share promoter 
elements. Previous work has shown that divergent gene- 
pairs are somewhat correlated in expression in yeast [21] 
and in several other species [22]. We have thus asked 
whether these pairs will also show correlation in mRNA 
degradation. Interestingly, we found statistically signifi- 
cant, yet modest, correlations between mRNA degrada- 
tion rates of divergent gene pairs, both with and without 
stress, and in the changes of degradation rates upon 
stress [10] (Pearson correlation = 0.24, 0.21, 0.25 and 
P-value = 2 x 10' , 1 x 10' , 1 x 10' , respectively, in Fig- 
ure S4 Additional file 1). We focus the following analysis 
in this section on the changes in mRNA degradation 
upon stress, as these reflect the condition-dependent reg- 
ulation of mRNA degradation, while absolute rates of 
mRNA degradation may be more dependent on intrinsic 
features of the mRNA (for example, secondary structure) 
that are not linked to transcription. Nonetheless, similar 
analysis of absolute mRNA degradation rates (in both 
conditions) revealed similar patterns (Figures S5 and S6 
in Additional file 1). 

Sharing of upstream intergenic regions among divergent 
gene pairs might not indicate shared promoters, especially 
when the size of the intergenic region is relatively large 
[23]. Moreover, the correlated regulatory patterns of 
neighboring genes could have multiple other causes, 
including chromosomal proximity, which might cause 
adjacent genes to be co-regulated by various mechanisms, 
non-random arrangement of microarray probes [24,25] 
and the evolutionary tendency of functionally related (and 
thus often co-regulated) genes to be clustered in the gen- 
ome [26]. However, these mechanisms should affect 
neighboring gene pairs of all orientations (Figure 3a), 
while promoter sharing should be more specific to diver- 
gent pairs. We have thus performed similar analysis for 
gene pairs of all three orientations to isolate the effect of 
promoter sharing from the additional factors that relate to 
genomic proximity (Figure 3b,c). Indeed, consistent with 
an effect of promoter sharing, gene-pairs with tandem and 
convergent architectures showed lower similarity in 
mRNA levels [21] and in mRNA degradation rates com- 
pared to divergent pairs. We further reasoned that promo- 
ter sharing is most likely in divergent gene-pairs that have 
short intergenic regions and high similarity in their mRNA 
expression profiles (and hence are likely co-regulated tran- 
scriptionally). We thus binned gene-pairs according to 
the length of their intergenic regions and their correlation 
in mRNA expression across a wide range of conditions 
[27] and examined the similarity in mRNA degradation 



profiles. As expected, the correlation in mRNA degrada- 
tion rates between divergent gene pairs increased with 
decreasing intergenic length and with increasing expres- 
sion similarity, and over all bins these correlations were 
more pronounced for divergent than for tandem or con- 
vergent gene pairs (Figure 3b,c). 

We also examined the similarity in mRNA degradation 
as deduced from the more recent and completely different 
technology of metabolic labeling [9], in which transcrip- 
tion is measured directly along with mRNA levels, and 
mRNA decay is deduced. This methodology was used to 
evaluate the changes in transcription and in mRNA degra- 
dation in response to an environmental stress. Interest- 
ingly, also in this dataset we detected significant similarity 
between the changes in mRNA degradation rates of diver- 
gently transcribed gene pairs, and a lower similarity of 
convergent and tandem pairs (Figure 3d,e). In addition, 
the same dependency on the length of the intergenic 
region and transcriptional similarity is also observed for 
divergent genes in this dataset (Figure 3d,e). Notably, the 
extent of correlated mRNA degradation among divergent 
gene pairs was only slightly lower than the extent of corre- 
lated transcription profiles (Figure S7 in Additional file 1). 

To add more support to the idea that shared promoter 
regulation results in coordination of post-transcriptional 
regulation of mRNA degradation, we examined pairs of 
genes that are not adjacent on the genome but share 
multiple promoter-associated regulators. For this we 
used physical measurements of protein-DNA interactions 
[28] and took all pairs of genes that share at least three 
regulators, and as a control we took pairs of genes that 
do not share a single regulator. We compared the differ- 
ential effect of oxidative stress on mRNA degradation 
between the transcriptionally co-regulated gene pairs and 
the control gene pairs, and found that transcriptionally 
co-regulated gene pairs have significantly more similar 
patterns of mRNA degradation compared to the control 
gene pairs (Figure S8 in Additional file 1). Note that we 
controlled for the expression similarity [10] of the gene 
pairs and therefore this result is not due to the possibility 
that genes with similar expression patterns are more 
likely to have similar patterns of mRNA degradation. 
Instead, the increased similarity in mRNA degradation 
appears to be linked to the co-regulation by promoter- 
associated transcription factors. 

In summary, although these analyses are correlative, 
they are based on patterns of hundreds of genes. Together 
they suggest that promoter sharing may contribute to the 
similarity in mRNA degradation profiles among neighbor- 
ing genes, and provide support to the notion that promo- 
ters regulate mRNA degradation. Note that correlated 
mRNA degradation is observed not only for divergent but 
also, more weakly, for convergent and tandem gene pairs. 
This might indicate that some regulatory elements could 
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regulate transcription and mRNA degradation from the 
terminator region, or alternatively that genomic proximity 
might enable sharing of transcription and degradation fac- 
tors that are recruited by regulatory elements of one gene 
and then co-regulate the neighboring gene. 



Coupling of transcription and mRNA degradation in 
mammals 

While the coordination of transcription and mRNA 
degradation has been studied primarily in yeast 
[6,7,11,12,17,18,29], previous studies also supported the 



Dori-Bachash et al. Genome Biology 2012, 13:R1 14 
http://genomebiology.com/201 2/1 3/1 2/R1 1 4 



Page 8 of 14 



possibility of a similar coordination in mammals. First, 
replacing the promoter of |3-globin was shown to affect 
its mRNA degradation [30]. Second, deletion of a tran- 
scription regulator (Matl of the TFIIH complex) led to 
compensatory changes in transcription and mRNA 
degradation [31]. Third, genes with a 'peaked' response 
to environmental changes were shown to have concomi- 
tant changes in transcription and mRNA degradation 
[13,14]. To further examine this possibility we analyzed 
the changes in mammalian mRNA degradation rates 
between different cell types or species. 

We first examined the differences in mRNA degrada- 
tion rates between induced pluripotent stem (iPS) cells 
and the human foreskin fibroblasts (HFFs) they were 
derived from, as recently measured by Neff et al. [32] . If 
changes in mRNA degradation upon induction of pluri- 
potency are independent of the transcriptional changes 
in this process, then we would expect to see a negative 
correlation between changes in mRNA degradation rates 
and in mRNA levels, since increased mRNA degradation 
leads to decreased mRNA levels. However, consistent 
with the possibility that, as in yeast, changes in mRNA 
degradation are coupled to larger changes in transcrip- 
tion, we observed a significant positive correlation 
between changes in mRNA degradation rates and in 
mRNA levels, such that increased iPS cell mRNA degra- 
dation is preferentially associated with increased (rather 
than decreased) iPS cell mRNA levels (Figure 4a). This 
result reproduces similar observations in yeast [10,12] 
and suggests a global coupling between human transcrip- 
tion and mRNA degradation that is manifested in the 
induction of pluripotency. 

We next performed two control analyses (Figure S9 in 
Additional file 1). First, in the above analysis the changes 
in mRNA levels and in mRNA degradation were calcu- 
lated from the same experiment and this may introduce 
spurious correlations. To avoid potential technical correla- 
tions, we instead compared the mRNA degradation 
changes as defined by Neff et al. to changes in mRNA 
levels as defined by the average of four other (indepen- 
dent) studies that each compared pluripotent to fibroblast 
human cells [33] . Second, the rate of decrease in mRNA 
levels following transcriptional arrest is assumed to reflect 
only mRNA degradation but could also involve other 
effects, such as residual transcription or differential degra- 
dation along the time course, which would confound the 
estimation of mRNA degradation rates. While we cannot 
completely exclude such confounding effects, we reasoned 
that these are less likely to have an impact on genes in 
which the decrease in mRNA levels upon transcriptional 
arrest has the smallest deviation from the expected model 
of exponential decay. We thus restricted the control analy- 
sis to 12% of the genes with the best fit to the expected 
mRNA degradation model. In both controls we observed a 



strong positive correlation between changes in mRNA 
degradation and in mRNA levels (R > 0.3, P < 0.001; 
Figure S9 in Additional file 1). 

Next, following our previous work in yeast [12], we 
sought to examine the possibility of transcription-degra- 
dation coupling from an evolutionary perspective 
through interspecies comparison of transcription and 
mRNA degradation. To this end, we compared mRNA 
degradation rates between human B cells (BL41) and 
murine NIH-3T3 fibroblasts [34]. The two datasets were 
produced by the same lab with the same methodology, 
but for non-orthologous cell types and therefore the dif- 
ferences reflect both evolutionary and developmental 
changes in gene regulation. This analysis is complemen- 
tary to the above analysis of HFFs versus iPS cells as it 
involves a different comparison of cell lines, but most 
importantly as it is based on a completely different 
methodology of metabolic labeling of newly transcribed 
mRNA with 4-thiouridine for 1 hour. Thus, transcrip- 
tion rates were estimated by the amount of newly tran- 
scribed mRNA, and degradation rates were estimated by 
comparing the newly transcribed to total mRNA levels 
(a high relative fraction of newly transcribed mRNAs 
indicates rapid turnover). 

We calculated interspecies differences in total mRNA 
levels (Atotal) and in newly transcribed mRNA (Anew), 
both reflecting log 2 -ratios of mouse versus human. Atotal 
and Anew were highly correlated (R = 0.88; Figure S10 in 
Additional file 1), indicating that most interspecies differ- 
ences in total mRNA levels are due to transcriptional differ- 
ences, or to differences in microarray probe efficiencies for 
orthologous human and mouse genes. Nonetheless, many 
genes had significant differences between Atotal and Anew, 
indicating differential mRNA degradation rates (for exam- 
ple, if Atotal > Anew, then the change in total mRNA levels 
cannot be accounted for by changes in transcription and 
suggests a change in mRNA degradation). To understand 
the connection between changes in mRNA degradation 
and changes in transcription, we considered three models, 
and the expected distributions of Atotal and Anew accord- 
ing to these models are shown for genes with increased 
(red) or decreased (green) mRNA degradation, in mouse 
compared to human (Figure 4b). Briefly, if differences in 
transcription are independent of the differences in mRNA 
degradation (Model I), then we expect no effect (equal dis- 
tribution for the two gene sets) in transcription (Anew), but 
a significant effect of mRNA degradation on total mRNA 
levels (Atotal). If differences in mRNA degradation primar- 
ily reflect measurement errors in Anew and Atotal (Model 
II), then we expect moderate shifts in the distributions of 
both Anew and Atotal (for example, increased mRNA 
degradation, Anew - Atotal > 0, would be identified prefer- 
entially at genes with Anew >0 and Atotal<G), as observed 
in comparison of unrelated (non-orthologous) genes 
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Figure 4 Human-mouse comparison suggests conservation of the coupling between transcription and mRNA degradation (a) Scatter 
plot of cell type-dependent differences in mRNA levels (y-axis) and in mRNA degradation rates (x-axis) for 4,335 genes, quantified as log 2 
(iPS/HFF). The Spearman correlation between the differences in mRNA degradation and mRNA levels is 044 and the corresponding linear least 
square fit is shown, (b) Theoretical distributions of Atotal (left) and Anew (right) according to three models, as detailed in the text, for all genes 
and for genes with differential mRNA degradation rates (red and green correspond to higher and lower degradation rates in mouse, 
respectively), (c) Observed distributions of Atotal (left) and Anew (right) following the color definitions in (b); vertical lines mark the distribution 
medians (Atotal: -0.16, 0.19 and 0; Anew: 1.33, 1.21 and 0; for red, green, and black, respectively). M.m and H.s denote Mus musculus and Homo 
sapiens. 
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(Figure Sll in Additional file 1). Finally, if, as in yeast 
[10,12], differences in mRNA degradation are typically 
coupled to larger differences in transcription (Model III), 
then we expect both Anew and Atotal to be shifted such 
that increased degradation is associated with increased 
newly transcribed and total mRNAs. 

Interestingly, comparison of human-mouse orthologs 
shows that the distributions of Atotal are, in fact, quite 
similar for genes with increased and decreased mRNA 
degradation, and are only slightly shifted (Figure 4c). As 
both of the first two models (Independent transcription or 
Noise) predict a large shift of Atotal in the same direction, 
the small shift in Atotal suggests compensation by tran- 
scriptional effects according to the third model (Coupled 
transcription), which predicts an opposite shift. Consistent 
with this possibility, we observe a remarkable shift in the 
distribution of Anew. The shift in Anew is much larger 
than the shift in Atotal, which can only be explained by 
the third model. These results indicate that, on average, 
differences in mRNA degradation are compensated for by 
transcriptional effects and may be explained by a combina- 
tion of the above models - in some genes transcriptional 
effects are independent of (and do not compensate) the 
changes in mRNA degradation, while in other genes tran- 
scriptional effects over-compensate for the changes in 
mRNA degradation and generate an opposite change of 
total mRNA levels. 

Discussion 

We previously predicted (using interspecific hybrids) cis- 
dependent evolutionary changes in mRNA levels and in 
mRNA degradation of hundreds of yeast genes. In the 
first part of this work we have validated these predictions 
and demonstrated that the relevant mutations often 
reside at the swapped upstream regulatory regions con- 
taining promoters and 5' UTRs. We further analyzed the 
contribution of different upstream segments to the diver- 
gence of two genes and identified promoter and UTR ele- 
ments that affect transcription and mRNA degradation. 
Promoter alterations could influence the recruitment of 
RNA-binding factors, such as Rpb4/7, that would shuttle 
with the mRNA to the cytoplasm and affect their degra- 
dation and/or translation. Other evidence for this model 
comes from our observation that gene pairs that share a 
promoter tend to have similar patterns of mRNA degra- 
dation across the yeast genome. This suggests that pro- 
moter-mediated effects on mRNA degradation often act 
in both directions of a bidirectional promoter. 

Recent studies identified two mechanisms by which 
promoter elements affect mRNA degradation [17,18] but, 
as noted above, these mechanisms are unlikely to account 
for our results. This may indicate the existence of multi- 
ple independent mechanisms for coupling transcription 
to mRNA degradation. Alternatively, it is possible that 



these seemingly distinct cases in fact reflect a common 
general coupling mechanism whose activity is modulated 
by multiple different regulators. Accordingly, general 
coupling factors (for example, Rpb4/7 and/or Ccr4-Not) 
may be recruited to a large fraction of the transcribed 
genes by yet unknown mechanisms and their gene- speci- 
fic activity could then be fine-tuned by various transcrip- 
tion regulators, including Rapl and cell-cycle regulators. 
This model is consistent with the various observations of 
a global (rather than gene-specific) coordination of tran- 
scription and mRNA degradation, as shown by this and 
previous studies [6,11,12]. This model is further sup- 
ported by observations that the same gene often displays 
a coupling of transcription and mRNA degradation 
through multiple distinct mechanisms. First, Bregman 
et al. [18] noted that apart from the large effect of Rapl 
binding sites, two other promoter elements (Pacti. 
Prpl3o) also had independent effects on the mRNA 
degradation rate of their construct. Second, our analysis 
of YLR326W identified distinct ds-regulatory elements 
that affected both transcription and mRNA degradation. 
Third, we found that genes in which transcription and 
mRNA degradation were previously shown to be coupled 
by cw-mutations were also coupled through the insertion 
of a selection marker. It would therefore be interesting to 
examine if Rpb4/7 and Ccr4-Not interact differently with 
mRNAs of coupled versus uncoupled genes. 

Finally, we analyzed differences in mRNA degradation 
between mammalian cell-types or species and provide 
evidence for a global coordination between transcription 
and mRNA degradation. Importantly, the analyses of 
mammalian cells, as well as the analyses of divergent 
yeast gene pairs, were both performed on two different 
datasets of mRNA degradation rates, with one that used 
transcriptional arrest and another that used metabolic 
labeling. These completely different technologies produce 
distinct estimates of mRNA degradation rates, yet in our 
analyses they show consistent results, further supporting 
our conclusions. Our observation of similar transcrip- 
tion-degradation coupling in yeast and in mammals may 
suggest that coupling is an evolutionarily conserved phe- 
nomenon. Given the conservation of the protein com- 
plexes Rpb4/7 and Ccr4-Not, it is tempting to speculate 
that similar mechanisms govern this coupling in yeast 
and in mammals and therefore that a promoter-mediated 
coupling exists also in mammals. Consistent with this 
possibility, recent work in a human embryonic cell line 
identified several mRNA-binding transcription regula- 
tors, as well as an overall enrichment of interactions 
between mRNA-binding proteins and transcriptional reg- 
ulators, demonstrating the many potential links between 
transcription and post-transcriptional regulation [35]. 
Alternatively, it is possible that the observed coupling 
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evolved independently in yeast and in mammals and 
therefore reflects convergent evolution. 

Conclusions 

Our results demonstrate, at the genomic scale, that tran- 
scription and mRNA degradation are coordinated in yeast 
as well as in humans. In yeast, detailed experimental and 
computational analysis further indicates that this coupling 
is driven by promoters and cannot be accounted for by 
two recently described mechanisms, suggesting that multi- 
ple mechanisms influence the coupling of transcription 
and mRNA degradation. Additional work would be 
needed to decipher the relevant mechanisms and their 
implications. 

Materials and methods 

Selection of genes for swapping experiments 

We first selected 40 genes for swapping experiments based 
on a high degree of inter-species as-dependent differences 
in mRNA levels and consistency of these cis-effects among 
several microarray experiments [36]. The first half of the 
genes (20 genes) was selected among the genes with a pre- 
vious observation of transcription-degradation coupling 
[12] and was selected to have equal representation of OPN 
and DPN promoters [37]. Of these 40 genes, six were dis- 
carded due to technical problems in generating the swap 
strains, and for the rest we examined the effect of swap- 
ping on mRNA levels. Of these, 10 genes with a previous 
observation of transcription-degradation coupling were 
selected at random and for these we examined the effect 
of swapping on mRNA degradation. One gene was 
excluded due to lack of reproducibility, and for the rest we 
estimated the swapping effect. Of these, the four genes in 
which the swapping effects on mRNA degradation 
appeared most reliable were chosen for further validation 
with a longer time course of eight time points. The swap- 
ping effect of one of these genes (GAL83) was largely 
reproduced in a control strain in which we only inserted 
the selection marker, while the swapping effect of another 
gene (CEX1) was relatively weak, and thus we further 
examined in detail (with a series of swap strains) the two 
remaining genes, MRU and YLR326W. The selection of 
swapped segments was based on the ability to design effi- 
cient primers, the location of neighboring genes and the 
location of known ds-regulatory motifs. 

Generating swap strains 

Wild- type S. cerevisiae BY4741 was compared to BY4741 
strains in which we replaced the promoter and 5' UTR 
segment of an individual gene with the orthologous seg- 
ment from S. paradoxus CBS432 (Table SI in Additional 
file 1). The orthologous segments were first cloned and 
inserted into a pBS7 plasmid next to a kanamycin selec- 
tion marker, and the resulting construct (selection marker 



upstream of a regulatory segment) was inserted into 

S. cerevisiae. 

In the first set of 17 genes, including all 9 genes for 
which we examined mRNA degradation, the S. cerevisiae 
upstream regions were replaced by this construct. For 
these genes we also generated control strains in which 
only the kanamycin marker was inserted (at the same posi- 
tion, but without replacing the regulatory region) and for 
two of these genes (YLR326W and MRU) we generated a 
series of strains where different portions of the upstream 
regulatory regions were replaced, thus generating chimeric 
intergenic regions composed of partial 5. paradoxus and 
partial S. cerevisiae segments (together encompassing the 
original length of the upstream region) and an upstream 
selection marker. 

A caveat of the above swapping approach is that the 
swapped region might contain an element (promoter or 
terminator) that regulates the adjacent upstream gene (not 
the one being studied), and these elements are not only 
swapped but are also being separated from the adjacent 
gene by the selection marker. To avoid this potential 
effect, in the second set of 17 genes we did not replace the 
S. cerevisiae regulatory region but instead inserted the 
construct downstream of the S. cerevisiae regulatory 
region (directly upstream of the start codon), thereby 
maintaining the S. cerevisiae regulatory region next to the 
adjacent gene but also inserting the S. paradoxus regula- 
tory region next to the gene being studied. For these genes 
we did not generate control strains. For measurements of 
mRNA degradation rates, the swap and control S. cerevi- 
siae strains were mated with a wild-type 5. paradoxus 
(CBS432) in order to allow comparison of the 5. cerevisiae 
allele (whose regulatory sequences are swapped) with the 
orthologous S. paradoxus allele, and also because tran- 
scriptional arrest experiments appeared to be more effi- 
cient in diploid strains (not shown). 

Comparison of wild-type and swap strains 

Real time PCR was performed for the WT and swap 
strains with duplicates or triplicates. We compared the 
mRNA levels of the genes whose regulatory sequences 
were swapped to reference genes that are not expected 
to be affected by the swapping and thus provide a nor- 
malization factor for comparison between the WT and 
the swap strains. The effect of swapping was thus quan- 
tified as: 

Aswap = log 2 [(SWg/WT g ) / (SW C /WT C )] 

where SW and WT denote the swapped and wild-type 
strains, respectively, and g and c denote the examined 
gene and its reference, respectively. ACT1 was used as a 
reference gene for PCR normalization in the analysis of 
mRNA levels (Figure lb). Preliminary analysis suggested 
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that ACT1 may be affected by the phenanthroline treat- 
ment and thus we did not use it as a reference in any of 
the mRNA degradation experiments. In the initial 
mRNA degradation experiments (Figure lc; Figure Sib 
in Additional file 1) we used ARO80 as a reference 
gene. We note that our experiments are designed to 
measure changes in mRNA degradation due to swapping 
of a single m-regulatory segment (rather than absolute 
mRNA degradation rates) and hence are less dependent 
on the choice of reference gene as long as the reference 
gene is not directly affected by the swapped regulatory 
segment. Nonetheless, for the validation experiments 
(Figure Id) we sought to use different reference genes in 
order to demonstrate the reproducibility of the observa- 
tions and their lack of dependence on the reference gene. 
For this reason and also in order to reduce the overall 
number of experiments, we used a circular design whereby 
each of the four genes being validated also served as a 
reference gene for one of the other genes: each gene was 
measured in the WT and in two swap strains including 
the swap of that gene and of another gene, where it served 
as a reference (GAL83 - YLR326W and CEX1 - MRU 
controlled for one another). 

Analysis of differential mRNA degradation 

Transcription was arrested by addition of 150 uM 1,10- 
phenanthroline to the WT and swap hybrid strains and 
monitored mRNA levels at four or eight time points. Least 
squares linear regression was used to estimate the slope 
along the time-course, which reflects the effect of swap- 
ping on rates of mRNA degradation. For example, a slope 
of -0.05 indicates that after 20 minutes the relative mRNA 
level of the swap strain decreases by 2-fold (log 2 [SW g / 
WT g ] = -0.05 x 20 = 0.5). These slopes were also used to 
estimate the relative effect of the swapping on mRNA 
degradation rates (shown in Figure 2) by comparison with 
the average half-life (ti/ 2 ) of the native S. cerevisiae gene 
from two previous studies [10,38]: 

A deg = log 2 [(l/ti/2 - slope) / (l/ti/ 2 )] . 

To determine significance of the differences in mRNA 
degradation rates between WT and Swapped (Figure 1), or 
between swapping of different lengths (Figure 2), we per- 
formed paired f-tests, evaluating the differences between 
the changes in mRNA levels at 20 minute intervals. In 
each experiment we can quantify the mRNA degradation 
rate at three 20 minute intervals: 

Ay+20 = log 2 (m i + 20 /mi) 

where mj is the relative mRNA level at time i compared 
to the reference gene, for i = 0,20,40. The £-test is per- 
formed over the differences between paired effects from 
the two strains (for example, Swapped and WT): A i>i+2 o 



(SW,k) - A u+20 (WT,k), for i = 0,20,40 and k = l..n, where 
n is the number of biological repeats. The pairing of corre- 
sponding time intervals is important because the pattern 
of changes along the time course is not exactly linear and 
is often gene-specific. By comparing the same gene 
between different strains and restricting the comparison to 
corresponding intervals we can largely eliminate the varia- 
bility along the time course and focus on the variability 
between strains (that is, the swapping effect). 

mRNA degradation experiments were typically per- 
formed with three technical replicates and with two biolo- 
gical replicates, with three exceptions: (1) YLR326W swaps 
of 200, 386, 398, 414, 444, 459 and 471 bp (Figure 2) were 
done in four replicates; (2) complete swaps of YLR326W, 
CEX1, GAL83 and MRU were analyzed without biological 
replicates but in two different experiments - an initial 
experiment with four time-points (Figure lc) and a later 
experiment with eight time points (Figure Id); (3) the com- 
plete swap of PEX32 was analyzed without biological repli- 
cates due to the strong effect of the selection marker. 

Human-mouse comparison 

Total and newly transcribed mRNA levels for human 
BL41 and murine NIH-3T3 cell lines with three replicates 
were taken from Gene Expression Omnibus (GSE10026). 
We averaged the three replicates, and compared the 
averages of selected probes for orthologous genes (follow- 
ing the orthologs and probe selection of Friedel et al. 
[34]). We calculated log 2 (mouse/human) for total (Atotal) 
and nascent (Anew) mRNA levels and centered these log 2 - 
ratios on zero. Human and mouse half-lives were taken 
from the supplement of Friedel et al. and genes with dif- 
ferential mRNA degradation rates were identified as those 
with at least 1.75-fold difference between the species both 
according to the Friedel et al. half-life estimates and 
according to (Anew - Atotal). We filtered out genes with 
very low hybridization intensities and/or very large inter- 
species differences, although the trends shown in Figure 4 
were maintained without filtering or with different strin- 
gencies of filtering (not shown). This analysis included 
4,209 genes (shown in black in Figure 4), with 729 genes 
defined as higher mRNA degradation rates in mouse (red) 
and 629 genes defined as higher mRNA degradation rates 
in human (green). 

5' RACE 

RNA was extracted from WT and two mutant strains 
using Yeast MasterPure kit (Epicenter Biotechnologies, 
Madison, WI, USA), including DNAse treatment to 
reduce contaminations of genomic DNA. To discriminate 
between mature capped mRNA and mRNA degradation 
products, RNA samples were dephosphorylated using 
FastAP thermosensitive Alkaline Phosphatase (Thermo 
Scientific, Waltham, MA, USA). RNA was then cleaned 
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using RNeasy MinElute Cleanup kit (Qiagen, Venlo , 
Netherlands ). 5' Caps were then removed using Tobacco 
acid pyrophosphatase (TAP) enzyme (Epicenter Bio- 
technologies) to generate 5' monophosphate for mature 
intact mRNAs and RNA was cleaned using RNeasy kit. A 
unique RNA adaptor (UCUUUCCCUACACGACGCU- 
CUUCCGAUCUGCGC) was ligated using T4 RNA ligase 
(New England BioLabs, Ipswich, Ma, USA) with incuba- 
tion of 1.5 hours at 22°C, and RNA was cleaned using 
RNeasy kit. First strand cDNA was generated from total 
amount of RNA ligation product, random hexamer pri- 
mers and M-MLV reverse transcriptase. cDNA (2 ul) was 
used as template for each following PCR reaction. To 
estimate the length of the 5' UTR a reverse primer 
located downstream of the start codon and a forward pri- 
mer on the adaptor (ATGATACGGCGACCACCGA- 
GATCTACACTCTTTCCCTACACGACGCTCTT) were 
used together and products were run on a 1.5% agarose 
gel. Reverse primer sequences: CGGTAGAAGTGCAAG- 
TAATGG for YLR326W and for ATAAATTAACGGCT 
GTGGGTC. Primers were validated using an additional 
forward primer located on the ATG of each gene, ATGT 
CAGGGTTCATTAAGAGC for YLR326W and ATGT 
CGTTGGAAGCCATCGTC for YPR118W. 

Additional material 
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